Intuitive understanding of non-gaussianity in ekpyrotic and cyclic models 
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It has been pointed out by several groups that ekpyrotic and cyclic models generate significant non- 
gaussianity. In this paper, we present a physically intuitive, semi-analytic estimate of the bispectrum. 
We show that, in all such models, there is an intrinsic contribution to the non-gaussianity parameter 
/jvi that is determined by the geometric mean of the equation of state Wefc during the ekpyrotic phase 
and Wc during the phase that curvature perturbations are generated and whose value is 0(100) or 
more times the intrinsic value predicted by simple slow-roll inflationary models, j™^''*"^"^ = 0(0.1). 
Other contributions to /atl, which we also estimate, can increase l/ivz,! but are unlikely to decrease 
it significantly, making non-gaussianity a useful test of these models. Furthermore, we discuss a 
predicted correlation between the non-gaussianity and scalar spectral index that sharpens the test. 



I. INTRODUCTION 

Observations of the cosmic microwave background and 
large scale structure have provided strong evidence for 
a nearly scale-invariant spectrum of primordial adiabatic 
density fluctuations with perhaps a slightly red tilt [1, 2]. 
The key cosmological question is: how was this spectrum 
generated? Two mechanisms are currently known, in- 
flation [3-6] and ekpyrosis [7-11]. Inflation is a period 
of accelerated expansion following the big bang, charac- 
terized by a large Hubble parameter H and an equation 
of state w —1. Ekpyrosis is a period of ultra-slow 
contraction preceding the big bang, characterized by a 
small H and w ^ 1. One way of distinguishing the two 
mechanisms is by measuring the spectrum of primordial 
gravitational waves [7, 8, 12] whose amplitude, propor- 
tional to , is exponentially different in the inflation- 
ary versus ekpyrotic/cyclic models. A second approach, 
recently emphasized by several groups analyzing differ- 
ent types of ekpyrotic models [13-18], is by measuring 
the non-gaussian contributions to the density fluctuation 
spectrum; e.g. the bispectrum, parameterized by the pa- 
rameter /tvl [19], as well as the correlation between f^L 
and the scalar spectral tilt ris [18]. 

The purpose of this paper is to provide a physically 
intuitive semi-analytic estimate of the contributions to 
the bispectrum that makes clear why measuring both the 
non-gaussianity and the scalar spectral tilt is a good ap- 
proach for testing ekpyrotic scenarios of various types, 
including the cyclic model. We show that ekpyrotic and 
cyclic models generically produce non-gaussianity of the 
"local" type [20] with a net value of Jnl that is at least an 
order of magnitude greater than the value predicted by 
simple slow-roll inflationary models (/atl ^1)- The pre- 
diction is timely because this value is large enough to be 
detectable in near-future cosmic microwave background 
studies. The current bound reported by the WMAP col- 
laboration, -9 < fNL < 111 and 0.94 < Ug < 0.99 at 2a 
[1] already intrudes on the ekpyrotic range. The uncer- 
tainty in /jv L will improve modestly with further WMAP 
data, and then dramatically (A/jv^ < 5) with the forth- 
coming Planck satellite mission and large scale structure 



studies. 

The reason why the non-gaussianity in ekpyrotic/cyclic 
models is generically several orders of magnitude greater 
than in inflationary models can be traced to the difference 
in the equation of state during the period that density 
fluctuations are generated. In standard versions of both 
models, the density perturbation spectra have their ori- 
gin in scalar fields which develop nearly scale invariant 
perturbations while evolving along an effective potential 
V{4)i). During an inflationary phase, the potential must 
be nearly constant in order to obtain Winf ~ —1 or, 
equivalently, Cinf = ^{1+Winf) ^ 1. This means that the 
inflaton is a nearly free field with nearly gaussian quan- 
tum fluctuations. The non-gaussian amplitude depends 
on the deviation of the potential from perfect flatness or, 
equivalently, on how close the slow-roll parameter emf 
and its time variation are to zero. This intuitive argu- 
ment is consistent with the quantitative expression ob- 
tained by Maldacena [21], for example. By contrast, a 
negative, exponentially steep potential is required to ob- 
tain Wek ^ 1 or, equivalently, Cefc ^ 1 in an ekpyrotic 
phase, which means that the scalar fields have significant 
nonlinear self-interactions whose magnitude depends on 
how large €ek is- Because the magnitude of eek is 0(100) 
or more times larger than ei„/, the scalar field contribu- 
tion to the non-gaussianity - which we will call the "in- 
trinsic" part " is significantly larger for ekpyrotic/cyclic 
models [10, 13-15, 17, 18]. 

This intuitive argument only refers to the intrinsic 
contribution to the non-gaussianity, but this is enough 
to argue why ekpyrotic/cyclic models generically pre- 
dict I/atl] to be several orders of magnitude greater than 
the value in inflationary models. Even if the additional 
effects discussed below add or subtract from [/atl] in 
the ekpyrotic/cyclic model, obtaining a value less than 
one (that is, in the inflationary range) would only occur 
through accidental cancellations of independent terms at 
the two- or three-decimal level, which is highly unnatural 
[22]. Conversely, it is possible to add features to infla- 
tionary models (such as curvatons, non-standard kinetic 
energy density, etc., with certain parameters) that en- 
hance the non-gaussianity beyond I/atl] = 1- However, 
these are unnecessary embellishments, and, when added. 
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can produce virtually arbitrary /nl of either sign. 

The predictions can be further refined by taking ac- 
count of how the scalar field fluctuations are transformed 
into curvature perturbations. In inflationary models, the 
scalar field fluctuations directly produce curvature per- 
turbations that are growing modes in an expanding uni- 
verse. Consequently, the tiny intrinsic non-gaussianity 
in the scalar fields discussed above translates directly to 
a tiny non-gaussianity in the curvature fluctuations (if 
no embellishments are added). In ekpyrotic/cyclic mod- 
els, the curvature perturbations produced directly by the 
scalar fields are decaying modes. Hence, the process in- 
volves scalar field fiuctuations first producing growing 
mode entropic perturbations during the ekpyrotic phase 
and, then, converting them to curvature perturbations 
just before the bounce to an expanding phase [10, 11]. 
It is notable that the equation of state during this con- 
version (or Ec) can be quite different from the equation 
of state during the ekpyrotic phase (or eefe). Since the 
entire curvature perturbation is produced by this con- 
version, even the intrinsic contribution is necessarily af- 
fected by Cc, as well as eefe. In fact, as we will show 
below, the intrinsic f^i^ turns out to be proportional to 
their geometric mean. Hence, the magnitude of /jvl in 
ekpyrotic/cyclic models can be considerably less if the 
conversion takes places in a kinetic energy dominated 
phase when = 3 <C Cefe, say, rather than in the ekpy- 
rotic phase [18]. This accounts for why the predictions 
for /at l were found to be significantly greater in the new 
ekpyrotic model [13-15] compared to the cyclic model 
[18], although both models predict values much greater 
than the inflationary case. 

The conversion mechanism is also important for deter- 
mining the sign of the intrinsic contribution to /nl- For 
example, the cyclic model discussed in Ref. [18] produces 
either sign for the intrinsic /nl, whereas the sign is typ- 
ically negative for the case considered in Ref. [13-15]. 
Furthermore, the conversion from entropic to curvature 
perturbations necessarily introduces an additional contri- 
bution to the non-gaussianity, which will be analyzed be- 
low. We will see that there are substantial regimes where 
this contribution is smaller in magnitude compared to the 
intrinsic contribution, but also substantial regimes where 
it is larger and can even reverse the sign. In the latter 
case, the net non-gaussianity tends to be so large that it 
is already ruled out by existing experiments. 

This paper is designed to translate the intuitive dis- 
cussion above into simple semi-analytic estimates, mak- 
ing clear those aspects that are highly model dependent 
and those that are generic. The conclusions themselves 
are not so original; for the most part, they appear in 
earlier papers focusing on particular examples [10, 13- 
15, 17, 18]. Our intent here is more modest: to provide 
simple expressions that clarify their origin and physi- 
cal interpretation, so that the signiflcance of forthcoming 
non-gaussianity tests can be better appreciated. The or- 
ganization of the paper is as follows: in the next section 
we set up our notation while reviewing the context of the 



calculations performed later in the paper. In section III 
we derive the intrinsic non-gaussianity present in the en- 
tropy perturbation in ekpyrotic models. This allows us to 
present an order-of-magnitude estimate of the resulting 
curvature perturbation in section IV, before proceeding 
to discuss various conversion mechanisms in more detail. 
Section V contains the conclusions that we draw from 
our analysis. 



II. SETUP 

In the ekpyrotic/cyclic model, the 4d effective theory 
describing the universe around the time of the big bang is 
given by gravity coupled to two minimally coupled scalar 
flelds with potentials. (This is the case in heterotic M- 
theory [23], for example.): 
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The scalar held and Friedmann equations are given by 



and 



+ 3H(j)i + Fi,^. = 



(2) 



(3) 



where H = a/a, a denotes the scale factor, and Vi^^^ = 
(dVi/dcpi) with no summation implied. In a contracting 
universe, a growing mode is given by the entropy pertur- 
bation, namely the relative fluctuation in the two fields, 
defined (at linear order) as follows [24] 



Ss = ((^1 5(t)2 - 02 S(l)i)/&, 

where we have defined 



(4) 



(5) 



The entropy perturbation is gauge-invariant and it rep- 
resents the perturbation orthogonal to the background 
scalar field trajectory, as shown in Fig. 1; see Ref. [25] 
for its definition to all orders. Its equation of motion, up 
to second order in field perturbations, is [25] 



6s + 3H5s + (Vss + 3^2^ 6s 



+ -{6s'"f + 
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(1) 



5-V,s-9-]i6s^'^r 
a a I 



(6) 



where we have neglected a non-local term that is \inim- 
portant in ekpyrotic models [18]. Note that this is a 
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closed equation for the entropy perturbation. Here de- 
notes a derivative of the potential along the background 
trajectory and the successive derivatives of the potential 
with respect to the entropy field are given by 
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(9) 



The angle 9 of the background trajectory is defined by 

[24] cos(6') = 0i/(T, sin(6l) = 02/o- and thus 9 = ~Vs/&. 
The parameter e and the equation of state parameter w 
are related to the background evolution via 



3/, ^ 



(10) 



and the ekpyrotic phase can be characterized by e 1. 



III. INTRINSIC NON-GAUSSIANITY 

The entropy perturbations are generated during the 
ekpyrotic phase, when the potentials are given by 

2 

Y,Vi{<ki) =-Vie-i^-'^' -V2e-i'''^^\ (11) 

i=l 

i.e. the potentials are negative and steep. Here we con- 
sider ci = ci(0i) > and C2 = 02(^2) to be slowly vary- 
ing in time, while Vx and V2 are positive constants. The 
ekpyrotic phase quickly flattens the universe, so that we 
can assume a flat Friedmann-Robertson- Walker (FRW) 
background with line element 5s^ — —6t^ + a^(t)5x^. If 
the Ci are exactly constant and moreover \ci\ ^ 1, then 
the Einstein-scalar equations admit the scaling solution 

a={-ty/^ 0, = ^ln(-^/^2^O, 7 = Ei 

This solution describes a very slowly contracting universe 
with e ^ 1. During the phase in which the entropic 
perturbations are generated, ^ = 0, which corresponds 
to a straight background scalar field trajectory. In this 
case, we define 



02 = 701, 

and, with this notation, Ci = 7C2 and 

I7C1C2I 



2(1+72)- 



(13) 



(14) 



Results are not very sensitive to 7 so we take 7 = 0{1) 
throughout. 



It is useful to recast the evolution in terms of the adia- 
batic and entropic variables a and s. Up to unimportant 
additive constants (that we will fix below), they can be 
defined by 



a = 



s = 



i>101 + 0202 



(15) 
(16) 



Then we can expand the potential up to third order as 
follows [26]: 

Fefe = -yoe^<^[l + es2 + |e3/V], (17) 

where K3 is of 0(1) for typical potentials (the case of 
exact exponentials corresponds to /ts = —4-^^2/3). The 
scaling solution (12) can be rewritten as 

a{t) = {-tf" a = -^ln (- ^e^*) s = 0. 

(18) 

The intrinsic non-gaussianity in the entropy perturbation 
is produced (hiring the ekpyrotic phase and can be deter- 
mined from the equation of motion for the entropy field 
(6), which reduces to 

6s + 3HSs + VssSs + ^VsssiSs^^^y = 0. (19) 



The last term is 



«3 r- 
sss — 2^ 



(20) 



and thus the solution, at long wavelengths and up to 
second order in field perturbations, is given by [15, 18] 



(5s(i) = 5s(i)(i) + 5s(2)(t) 

c- ^end . ~( c ^end \2 
= OSend—^ + C{dSend^^) , 



where 



K3^/e 



(21) 



(22) 



Thus, the intrinsic non-gaussianity present in ekpyrotic 
models (second term in (21) above) is of order 0{^/e^) 
and it is due to the steepness of the potentials and the 
resulting self-interactions of the scalar fields. This is in 
sharp contrast with inflation, where the intrinsic non- 
gaussianity is extremely small due to the flatness of the 
potential (or, equivalently, ei„f <C 1). 



IV. CONVERSION 
A. An estimate of intrinsic /jvz. 

What is measured is not directly the non-gaussianity 

present in the entropy perturbation, but the non- 
gaussianity it imprints on the curvature perturbation. 
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Thus, it is important to know the strength with which 
this intrinsic non-gaussianity gets transferred to the cur- 
vature perturbation. The time evolution of the curvature 
perturbation to second order in field perturbations and 
at long wavelengths is given in FRW time by [25] 

IM FT V 

n = —ess+^[-{v,s+^e^){Ss^'^f + ^6sSs]. (23) 

a a a 

where we have omitted a non-local term which can be 
neglected in ekpyrotic models for the same reason as the 
non-local term omitted from (6) above, see [18] for de- 
tails. At linear order, a non-zero entropy perturbation 
combined with a bending (9 ^ 0) of the background tra- 
jectory sources the curvature perturbation on large scales 
and results in a linear, gaussian curvature perturbation 




(24) 
(25) 



where we denote the duration of the conversion by At and 
the sign corresponds to a contracting universe. Thus the 
strength of conversion is proportional to 1/ ^/e^ and this 
dependence on the equation of state will have repercus- 
sions for the magnitude of the second order correction as 
well. 

Here, as with all the contributions to non-gaussianity 
discussed in this paper, the fluctuations arc generated by 
the scalar fields with canonical kinetic energy density, so 
they generate non-gaussianity of the "local" type, as de- 
fined in Rcfs. [19-21]. The local wavelength-independent 
non-gaussian contribution to TZ can then be character- 
ized in terms of the leading linear, gaussian curvature 
perturbation TZl according to 
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(26) 



using the sign convention for wavelength-independent 
non-gaussianity parameter /jvl in [19]. Then the con- 
tributions to /nl can be divided into three parts [18] 



NL 



'•^^L J At 



f-SsSs] 

a 

rintegrated _ __^_ls:Jl) U ,\\2 
JNL ~ Vl'R? y^end)) 



(27) 

(28) 
(29) 



jvntrmstc ^riscs from the direct translation of the intrin- 
sic non-linearity present in the entropy perturbation into 
a corresponding non-linearity in the curvature perturba- 
tion. By contrast, f^f^^t^"" and f^l^arated ^^^j^ 
non-zero even if the entropy perturbation were exactly 
gaussian. Both contributions are due to the non-linear 



relationship between the curvature perturbation and the 
entropy perturbation, as expressed in equation (23) - the 



difference is that f^^j^ 



integrated 



gets generated during the 



ekpyrotic phase and J^^^^*^"^ during the conversion. 

The intrinsic non-gaussianity can be estimated by com- 
bining equations (21), (25) and (27): 



^intrinsic 
JNL 



± 



3111 



At V ec 



(30) 



It is given by the geometric mean of the e parameters 
during the phases of generation and conversion of the 
perturbations. This is perhaps the most important result 
because it fixes a rough magnitude for IJnl] which can 
only be significantly reduced through accidental cancella- 
tions due to the other terms or other effects not included 
in the present analysis. It explains in a nutshell why the 
non-gaussianity in ekpyrotic/cyclic models is necessarily 
more than an order of magnitude greater than in sim- 
ple inflation models, and it explains why the equation of 
state during conversion can have a significant quantita- 
tive effect on the prediction. 

To go beyond this qualitative estimate to a more pre- 
cise one, we need to take account of the details of the 
conversion mechanism. We will now discuss the various 
possibilities considered in the literature, namely conver- 
sion during the ekpyrotic phase, after the ekpyrotic phase 
during kinetic energy domination, and conversion after 
the big bang by modulated preheating. 



B. Conversion during kinetic energy domination 

In the original ekpyrotic and cyclic models, the phase 
dominated by the steep, ekpyrotic potential V{(j)) comes 
to an end (at t = tend < 0) before the big crunch/big 
bang transition (i = 0), and the universe becomes dom- 
inated by the kinetic energy of the scalar fields. Conse- 
quently, the equation of state a.t t = tend changes from 
EeA: ^ 1 to e = 3 (corresponding to w; ^ 1, the equation 
of state for a kinetic energy dominated universe). In this 
subsection, we consider the case where the conversion 
from entropic to curvature perturbations occurs during 
this kinetic energy dominated phase. This occurs natu- 
rally in the heterotic M-theory embedding of the cyclic 
model because the negative-tension brane bounces off a 
spacetime singularity [27] - creating a bend in the tra- 
jectory in field space in the 4d effective theory - before it 
collides with the positive-tension brane (the big crunch- 
big bang transition). This bending of the trajectory auto- 
matically induces the conversion of entropy to curvature 
perturbations [11]. Here we will not restrict the analysis 
to this particular example, but we will consider the gen- 
eral situation in which the scalar field trajectory bends 
in a smooth way during the phase of kinetic energy dom- 
ination following an ekpyrotic phase. Before proceeding, 
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FIG. 1: The trajectory in field space reflects off a boundary 
at (f)2 = 0. The entropy perturbation, denoted 5s, is orthog- 
onal to the trajectory. The bending causes the conversion of 
entropy modes into adiabatic modes 5a, which are perturba- 
tions tangential to the trajectory. 



we should remark that due to the instability of the scal- 
ing solution during the ekpyrotic phase [11, 28], we must 
assume the background trajectory is localized very close 
to the ridge of the two-field ekpyrotic potential at the 
beginning of the ekpyrotic phase in order for the conver- 
sion not to happen already during the ekpyrotic phase. 
How this initial condition can arise in a cyclic model is 
beyond the scope of the present paper. 

We can immediately perform an order-of-magnitude 
estimate of the non-gaussianity since during the kinetic 
phase tc — 3, and so (30) would lead us to expect Jnl 
to be of order 



(31) 



We will treat the bend as if only one field reflects (02 )> 
and we will consider the case where > 0, see Fig. 1. 
Other cases can be related to these representative ex- 
amples by changing the coordinate system in field space 
appropriately. In the heterotic M-theory example, the 
reflection occurs because (j)2 comes close to a boundary 
of moduH space ((/)2 — 0) and is forced to bounce [27]. 
For the purposes of this study, we will treat the reflec- 
tion as being due to a potential, F^((/)2), an additional 
contribution unrelated to the exponential potentials that 
were dominant during the ekpyrotic phase (but which are 
negligible during the kinetic energy dominated phase). 

It is important to know the evolution of the entropy 
perturbation during the process of conversion. If there 
is a phase of pure kinetic energy domination before the 
conversion, then the background scalar field trajectory is 
also a straight line during this phase, but with the po- 
tentials being irrelevant. The equation of motion reduces 
to 



5s + i(5s = 0, 



(32) 



5s{t) = 5s(^\t) + 5s^^\t) 



= <5se„d(l + In ^) + c<5sL,(l + 2 In :^133) 



Incidentally, note that 



In(-i) 



(34) 



during the kinetic phase, while during the ekpyrotic 
phase tSs ~ 5s. This observation will simplify our analy- 
sis later on. 

During the conversion, even though the kinetic energy 
of the scalar fields is still the dominant contribution to 
the total energy, the potential V^'{(j)2) that causes the 
bending has a significant influence on the evolution of the 
entropy perturbation. To analyze this, we will approxi- 
mate (as in [16]) the bending of the trajectory to be grad- 
ual by taking 9 constant and non-zero for a period of time 
At, starting from t = tref- Note that, assuming the total 
angle of bending is 0(1) radian, we have |^| « l/|ire/| 
in this case. Then one can relate the derivatives of the 
potential to expressions involving 9, for example 



It' 



(35) 



Assuming a gradual conversion (At ~ ^re/), we can ig- 
nore higher derivatives of 9. (In [18] it was shown that 
sharp transitions lead to unacceptably large values of 
/atl; hence these cases are of less phenomenological in- 
terest.) To satisfy the constraint on the amplitude of the 
curvature perturbation obtained by the WMAP observa- 
tions, we set tref ~ —lO'^Mp^ [11]. At linear order, the 
equation of motion (6) then reads 



(1) 



3H5. 



(1) 



1 



9)5s 



(1) 



0, 



(1) 



(36) 



As indicated by equation (34), we can set 6s ' = as a 
first approximation and, thus, neglect the damping term 
in the equation of motion. Also, we will simply evalu- 
ate the coefficient of the last term midway through the 
reflection, and define 




1 



9j{tref + At/2) 



(37) 



For a gradual reflection a; « (2 — 3)^. Then the solution 
for the linear entropy perturbation is 



Ss^^"^ = 5s'--^\tref) COSUj{t — tref)- 



(38) 



and by matching onto the ekpyrotic solution and its first 



Instead of continuing to grow logarithmically, the entropy 
perturbation actually falls off during the conversion, see 
Fig. 2. This has the consequence that the conversion 
from entropy to curvature perturbations is less efficient 
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FIG. 2: The evolution of the Unear entropy perturbation dur- 
ing conversion: the sohd (blue) line shows the actual evolution 
calculated numerically, while the dashed (purple) line shows 
the approximate solution (38) with to = 6 = — 1/tre/ and 



re/ 



-400Mp/. For the purposes of illustration, 5s^^\t 



ref ) 



has been normalized to 1. 



than one might have naively thought. From (25), we can 
estimate 



7^L = 



.(1) - 



At 



3 uj 



ref 



(39) 

where we have used ec ~ 3, which is a good approxima- 
tion for subdominant reflections (and acceptable for the 
estimating purposes for dominant ones). Since ujAt « 3, 
the entropy perturbation evolves over nearly a half-cycle 
and consequently sinwAi is a small factor which we will 
take to be about 1/3 in our estimates, which fits well 
with numerical results. 

We also need to know the evolution of the second or- 
der entropy perturbation during the time of conversion. 
Using the same approximations as above (implying that 
we can neglect Vsss compared to 9Vss/ct and O^/ahy the 
time the conversion is underway) , the equation of motion 
(6) simplifies to 



(2) 



2^5(2) _V(5sW)2 = 0, (40) 

a 



where Ss'^^^ is given in (38). Putting ^ « tr at the start 
of the reflection and keeping in mind that we impose the 

• (2) 

boundary condition 5s ~ 0, the solution for the second 
order entropy perturbation is given by 

Js(2) ^ Ss'^^'>{tref)cOs[u{t-tref)] 

+ ^('5s^'HW))'( - 4cos[^(< - Uef)] 
-|-4cos'*[a;(< — tref)] + 9 sin^[cj(i — i^e/)] 
-I- sin[a;(< — tref)] sm[3uj{t — Uef)] ) . (41) 



At large eek (or, equivalently, large |c|), the second or- 
der entropy perturbation falls off in the same way as the 



linear perturbation, but at small |c| there are significant 
corrections to this behavior. We will also need the inte- 
gral 



c- 



^^(2) _ (1 + 2ln{tend/tref)) Sm{ujAt) ^ 
At (1 + lll(^end/W))^ ^ 

At sm{ujAt) sin(2tjAt) 

+ _ . 



(42) 



3uj 12lo 

where C = 5s'^^\tref)- In all cases of interest, the last 
two terms are negligible. 

We are finally in a position to evaluate the various 
contributions to the non-linearity parameter Jml- The 
intrinsic contribution, defined in (27), becomes 



intrinsic 
NL 



(43) 



with 



A 
B 



5UJ{1 + 2lll{tend/tref)) 

SVeeil + In (tend/tref))^ sm{ujAt) 
2VG0^ sm^{u;At)' 



(44) 



Eq. (43) is one of our key results because it shows that 
the essential contribution scales in a simple way with eek 
and has a value that exceeds the total f^L for simple 
infiationary models by more than an order of magnitude. 
As suggested by Eq. (30), the first term in J™2'"™szc ^^^^ 
be re-expressed as O (y/ec^ek) with ec = 3. This contri- 
bution comes directly from the non-zero 5s generated 
during the ekpyrotic phase, which is due to the third 
derivative w.r.t. s of the ekpyrotic potential. It is in- 
teresting to note that, in the case where Vsss = or 
K3 = 0, is nevertheless non-negligible because 

of the positive offset B generated by the linear entropy 
perturbation Js'^^ during conversion, the piece propor- 
tional to ((5s(i))2 in (41). 

Using the definition (28) together with (35) it is 
straightforward to see that for conversion during the ki- 

1 j-re flection • i i -i, u 

netic phase, /jy£ is always negative, and it can be 

estimated as 



pre flection 
JNL 



67^i J At 



{29h 



15a;" 



--)(,5s«)2 

7 

\tref + At/2| 



8^2 sin^icuAt) 

where we have used ^trefO ~ 1 and 

1 



At 



At. 



(45) 



tsin' ijit- tref) ^ -(.tref + ^) At, (46) 
At II 

which is a valid approximation since Js*-^^ evolves over 
approximately a half-cycle. We have neglected the con- 
tribution due to the term proportional to bsbs, since it is 
suppressed on account of (34). The integrated contribu- 
tion to /atl generated during the ekpyrotic phase gives 
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an additional contribution of 



/integrated 
NL ^ 



l27^i 



^ ( 1 + In {tend/tref W Sm (ujAt) 



(47) 



Note that neither the reflected nor the integrated contri- 
butions depend on eek', they both simply shift the final 
result by a number depending on the sharpness of the 
transition and the duration of the purely kinetic phase 
respectively. 

The total Jnl is the sum of all the above contributions. 
Since we are only considering gradual conversions, we ex- 
pect the ratio \tref + At/2\/At to lie between 1 and 2. 
Also, the duration of the pure kinetic phase between the 
end of the ekpyrotic phase and the conversion is necessar- 
ily rather short, so that we expect ln(tend/tref) ^ C'(l)- 

Putting everything together, we get the following esti- 
mates: 



f intrinsic 

Jnl 

pref lection 

Jnl 

pintegrated 
JNL 



IQc + 50 

-50 

5. 



(48) 
(49) 
(50) 



The analytic estimate for suggests a range that 

could extend to —100, but the value above is more con- 
sistent with the exact numerical calculations for a wide 
range of parameters. Altogether, we end up with the 
following fitting formula 



rtotal 
JNL 




(51) 
(52) 



Fig. 3 shows that this formula agrees well with the re- 
sults from exact numerical calculations. In particular, 
we have been able to estimate the slope, and hence the 
dependence on Ks^eefc, rather accurately. 



C. Conversion during the ekpyrotic phase 

In the "new ekpyrotic" scenario [13, 14, 29], the conver- 
sion of entropy to curvature perturbations takes place be- 
cause the background trajectory switches from the two- 
field unstable scaling solution to a single-field attractor 
solution; in other words, the trajectory starts out close 
to the ridge of the two-field potential and then falls off 
one of the steep sides, either because of the initial con- 
ditions or due to an additional feature in the potential. 
Adding a feature to the potential makes little difference 
to the results as long as the reflection remains gradual, 
so in fact we will not consider an additional potential 
in this section. Since the conversion happens during 
the ekpyrotic phase, i.e. while the ekpyrotic potentials 
are relevant, Cc = Cek, and from (30) we expect the non- 
gaussianity to be of order 




FIG. 3: A comparison of the results from numerical calcu- 
lations with the fitting formula given in Eq. (52) and indi- 
cated by the thick black line. Here we have fixed the value of 
Eefc = 36. The plot confirms that Jml then grows linearly with 
K3, in good agreement with (52). The sample models (dashed 
and dotted lines) are representative of the range of models 
and parameters shown in [18]. We have similarly checked the 
dependence on y/e when the value of K3 is kept fixed. 



As shown by Koyama et a/. [15], the 5N formalism is well 
suited to treating this case, with the result that the total 
local non-linearity parameter is given by 



NL 



12"^^' 



(54) 



where the index j corresponds to the field that becomes 
frozen in the late-time single-field solution (here we have 
reverted to the potential form (11), as it is better suited 
for studying conversion during the ekpyrotic phase). As 
a check, we have performed the direct integration of the 
equations of motion numerically [30] , and we find results 
that are in good agreement with the above analytic es- 
timate. Note that the sign of J^l is always negative for 
these cases. 

A qualitative understanding of this result can be 
achieved in the present context as follows: from (27) we 
can see that, since H/a is negative and approximately 
constant, the sign of J™^''™'**^ given by the sign of 

J 9Ss^^\ In practice, numerical simulations indicate that, 
for conversion during the ekpyrotic phase, we have 

^ (55) 



' ' 10|i| 

during most of the conversion, but with 9 growing faster 
towards the end. We will use this numerical input to 
guide our analysis. The sign of y™^''*"'**^ is essentially 
determined by the sign of 9Ss'-^'> towards the end of the 
period of conversion. Naively it is difficult to perform a 
purely analytic estimate of /nl in the current scheme, 
since all the terms in the equation of motion (6) go as 
Cii~'*; however, equation (55) tells us that 



fNL ^ eek ^ 0{cl). 



(53) 



(56) 
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Cl 


C2 


fNL,SN 


Jnl 


10 


10 


-41.67 


-39.95 


10 


15 


-41.67 


-40.45 


10 


20 


-41.67 


-40.62 


15 


10 


-93.75 


-91.01 


15 


15 


-93.75 


-92.11 


15 


20 


-93.75 


-92.49 


20 


10 


-166.7 


-162.5 


20 


15 


-166.7 


-164.4 


20 


20 


-166.7 


-165.1 



TABLE I: Ekpyrotic conversion: the values of /jvz, estimated 
by the ^A'^ formalism compared to the numerical results ob- 
tained by directly integrating the equations of motion. 



and so there are surprisingly few terms in the equation 

of motion for Ss^"^^ that arc actually important during 
(most of) the time of conversion. In fact, to a first ap- 
proximation, we are simply left with 



z a 



(57) 



Vsss decreases in importance as the single-field scaling 

solution is reached. Thus, even though the Vsss term 
determines the initial evolution of 5s^'^\ eventually the 
Vss term dominates. Then, since Vss < it is easy to 
sec that the sign of (^.s^^^ is always driven to be opposite 
to that of 9, and consequently y™^''"'**'^ is negative in all 
cases. 

Since j^<'/''='=**°" jg proportional to 7?.£^, it provides 
a contribution of the same order 0{eek)- Using (56), 
it is straightforward to see that the part of y^^/'^'^*'"" 
that is proportional to {5sY is always positive, while the 
part proportional to 5s6s is always negative. This im- 
plies that there will be a competition between the two 
contributions. Numerical integration then shows that 
the part proportional to ((5s) ^ approximately cancels out 
pntrmsic^ while the 6s5s part by itself is very close in nu- 
merical value to the final answer. In all cases f^^'^^rate 
is completely negligible. 

Consequently, the total /jvl turns out to be negative 
and moreover in good agreement with the SN result (54) 
[15], as shown in Table I. Clearly, the SN formalism 
provides a fast and elegant derivation of non-gaussianity 
for conversion during the ekpyrotic phase that agrees well 
with direct integration of the equations of motion. (Note, 
however, no analogous SN approac;li applies to conversion 
in the kinetic energy dominated phase.) 

The conclusion that /nl is substantially less than zero 
differs from Ref. [16]; their approximation method was 
incomplete though, since the term proportional to Ss5s 
in (23) was not included. As shown above, this term typ- 
ically contributes significantly to the final result. Also, 
the evolution of the entropy perturbation during conver- 
sion was neglected, which meant, for example, that val- 
ues of f'^^'trmstc qJ^^j^qj- gjgj^ were obtained. Our results 
here show that their approximations were too crude, and 



that Jml is generally negative with J^l ^ ^20. This is 
inconsistent with current limits by roughly 3 a [1]. 

Hence, we conclude that models with conversion during 
the ekpyrotic phase (tc = Cefe ^ 1 j e difficult to accom- 
modate with current observations in contrast to models in 
which conversion occurs during the kinetic energy domi- 
nated phase (cc = 3^. 



D. Conversion after the crunch/bang transition 

It has recently been proposed by Battefeld [17] that, 
instead of converting entropy perturbations into curva- 
ture perturbations before the big crunch/big bang transi- 
tion, the conversion could occur during the phase shortly 
following the bang through modulated reheating. The 
concept is that massive matter fields are produced co- 
piously at the brane collision and dominate the energy 
density immediately after the bang. The massive fields 
are assumed to couple to ordinary matter with a strength 
proportional to h(Ss), so that their decay into ordinary 
matter occurs at slightly different times depending on the 
value of ds. In this way, the ordinary matter perturba- 
tions inherit the entropic perturbation spectrum. Since 
the conversion happens while Ec ~ 3, we can expect an 
intrinsic contribution to Jnl of 



fNL ~ \/eIfc ~ ^(ci) 



(58) 



i-e. fNL is of the same order as in the case of conver- 
sion during the kinetic phase preceding the big crunch. 
Moreover, as shown in [17], no significant additional con- 
tributions to /jvL are expected. Therefore, models of 
this type are subject to roughly the same observational 
constraints as the models where the conversion happens 
during the kinetic phase before the bang and where the 
intrinsic contribution to the non-gaussianity is dominant. 

Note however that a detailed prediction is made dif- 
ficult by the fact that is an unknown function of 5s. 
Indeed, the entropy perturbations are converted with an 
efficiency [17] 



e = 



3 |/i 
2 



h 



(59) 



Since the non-linearity in the entropy field is of magni- 
tude c, this implies a non-linearity in the curvature per- 
turbation given by 



e 



(60) 



where the ± sign reficcts our ignorance of the sign of h^s 
and of the direction of bending of the scalar field trajec- 
tory. Thus, in the absence of a more detailed model, we 
cannot go beyond this rough order-of-magnitude estimate 
at present. 
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V. CONCLUSIONS 

Our analysis can be summarized in a few rules of 
thumb: 

• The intrinsic contribution to I/atlI is proportional 
to the geometric mean of Cek and ec, which is 
at least two orders of magnitude greater than in 
simple inflationary models (where | Z™!'^™'**'^ I = 
C(O.l)). When all contributions are considered, the 
total Jnl is generically more than an order of mag- 
nitude greater than in simple inflationary models, 
(where ' = 0{l)). 

• For a fixed value of K3, the value of is corre- 
lated with the spectral tilt: smaller implies 
smaller tek^ which tends to make the spectral tilt 
bluer. Current limits on f^L fit well with limits 
on the spectral tilt for the simplest models with 
conversion during the kinetic energy driven phase 
before the bang or during reheating after the big 
bang. 

• Models in which the conversion occurs during a 
phase with larger ec produce a larger intrinsic I/atlI 
and are more difficult to fit with the current ob- 
served limits on /^vi and spectral tilt. In particu- 
lar, models with conversion in the ekpyrotic phase 
favor |/jvl| large and negative. 

• Cases in which the intrinsic contribution to Jnl 
is much smaller than the reflection plus integrated 
contributions produces very large values of I/atlI ^ 
100 that are inconsistent with current observational 
bounds. This includes all cases where the conver- 
sion is sharp. 

The analysis suggests a useful characteristic plot for 
differentiating cosmological models: /atl versus tilt. Fig- 
ure 4 illustrates the prediction for a cyclic model in which 
the conversion from entropic to curvature perturbations 
occurs in the kinetic energy dominated epoch following 
the ekpyrotic phase. Here we keep the skewness K3 of the 
potential fixed, as we vary its steepness e. The prediction 
is a swath whose width is largely due to the uncertainty in 
the trajectory, parameterized by 7. Although the swath 
includes /jvi near zero, positive /jvl between 10 and 100 
is preferred for tilts in the range suggested by WMAP5. 
The prediction for simple inflationary models is confined 
to the narrow band A/jvl ^ 1 around zero. 

The results are surprisingly predictive. If observations 
of Jnl lie in the range predicted by the intrinsic contribu- 
tion of either inflationary or ekpyrotic/cyclic models, it is 
reasonable to apply Occam's razor and Bayesian analysis 
to favor one cosmological model over the other. Combin- 
ing with measurements of the spectral tilt significantly 



sharpens the test. The current observational bounds ob- 
tained by the WMAP satellite are still inconclusive, but 
it is clear from the estimates presented here that non- 
gaussianity should be detected by the Planck satellite if 




FIG. 4; A plot for characterizing the correlation between 
j/jv^l and scalar spectral tilt, — 1, here illustrated for the 
case of the cyclic model in which the conversion from entropic 
to curvature perturbations occurs during the kinetic energy 
dominated phase just before the big crunch/big bang transi- 
tion. Different curves correspond to different fixed amounts 
of skewness K3 in the potential (the central curve corresponds 
to = 4-\/2/3), while we vary the steepness of the potential 
Eefc. The curves show the general trend that j/jvi] increases 
as the spectrum becomes redder. Simple infiationary models 
correspond to the narrow horizontal hashed (red) strip with 
l/iVL| ^ 1. The shaded rectangle represents the current ob- 
servational constraints on /jvl and tilt (95% confidence) from 
WMAP5 [1]. 



the ekpyrotic/cyclic model is correct. At the same time, 
this provides a strong incentive to further refine other 
methods of measuring non-gaussianity, such as looking 
for evidence in measurements of the large scale structure 
of the universe. 
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